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NONLINEAR FLUID OSCILLATIONS IN A PARTIALLY FILLED 
AX I SYMMETRIC CONTAINER OF GENERAL SHAPE 

INTRODUCTION 

The motion of a fluid, having a free surface, has been the subject of 
much interest to men of inquiring minds throughout history. Problems such 
as the wave motion over sloping beaches, flood waves in the rivers, and the 
motion of ships in the sea preoccupied scientists in earlier days. Lagrange 
investigated surface wave problems and was a forerunner in this field to such 
noted scientists as Cauchy, Poisson, Airy, Stokes, Kelvin, Rayleigh, Lamb, 
and Stoker. These men made a great number and variety of studies of the 
physical problems involving surface waves and were outstanding in their 
contributions to this field. 

Research in the field of surface wave problems, in the last two 
decades, has produced more information than all preceding studies. One 
reason for the rapid pace of surface wave research in this period can be 
attributed to the advent of high speed computers which permitted the use of 
methods of analysis that were not practical earlier. A second reason is that 
the aerospace industry needed vast amounts of surface wave information to 
design large jet aircraft and space vehicles. 


This investigation is concerned with the type of surface wave problems 
known as standing gravity waves of finite -amplitude. W. G. Penny and A. I. 
Price [ 1] analyzed such waves in a rectangular tank of infinite depth and 
published the first study of finite -amplitude standing waves in 1952. Since 
then, a number of publications have appeared generalizing and extending their 
work. 

Tadjbakhsh and Keller [2] analyzed finite -amplitude standing waves 
in a rectangular tank, as did Penny and Price ; however, their analysis was 
for a finite depth tank. They found that finite standing gravity waves in a 
rectangular container will have a lower frequency than infinitesimal standing 
waves at large depths, but they have higher frequency below a certain depth. 

Fultz [3] , disturbed by the frequency effect reversal found by 
Tadjbakhsh and Keller, undertook an experimental study of finite -amplitude 
standing gravity waves in rectangular tanks. He found experimentally, as did 
Tadjbakhsh and Keller theoretically, that there is a frequency effect reversal 
of a fluid in a rectangular tank. 

Mack [4] analyzed standing gravity waves of finite -amplitude in a 
circular cylinder and obtained results for the first axisymmetric standing 
wave. DiMaggio and Rehm [5] then studied the same circular cylinder to 
obtain the first nonaxisymmetric standing wave . 
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Baird [6] investigated the first nonaxisymmetric standing wave in 
a cylindrical sector container. 

Moiseyev [7] published a general theory for free and forced oscilla- 
tions of a fluid in a rigid container. To the author's knowledge, this method 
has not been used to obtain finite -amplitude gravity waves because of the 
inherent difficulty of solving the integral equations upon which the theory is 
based. 

The problem set forth here is to obtain the first nonaxisymmetric 
finite -amplitude gravity wave in an axisymmetric rigid container. It is 
assumed that the fluid is perfect and surface tension effects are negligible. 

The formulation of the problem results in a nonlinear boundary value 
problem. The nonlinearity is in the boundary condition at the free surface 
of the fluid. Not only is the boundary condition at the free surface nonlinear, 
but also it must be satisfied on a moving surface of unknown shape. The 
method of solution is to expand the nonlinear boundary in a Taylor's series 
about the undisturbed free surface, and then satisfy it asymptotically by the 
method of Krylov and Bogoliubov [8] . 

The solution to this problem is of dual interest. Not only is it a 
fundamental study in the field of surface waves, but also it provides a means 
to obtain a thorough understanding of surface waves in space vehicle propel- 
lant tanks . In no other application to axisymmetric containers is the motion 
of the fluid more important than in the field of space research since the 
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propellant in a space vehicle contributes over 90 percent of the total weight. 
Thus, this study contributes much information about the overall dynamic 
characteristics of a space vehicle . 


4 



GOVERNING EQUATIONS FOR FREE OSCILLATIONS OF A FLUID IN 
A RIGID, PARTIALLY FILLED AX I SYMMETRIC CONTAINER 

The volume of fluid, V , in the axisymmetric container is bounded 
by the free surface, S , and the surface in contact with the container, F . 

The free surface existing when the fluid is at rest is known as the undisturbed 
free surface, S 0 . 

In this analysis the fluid is assumed to be incompressible and inviscid 
with its motion irrotational . The effects of surface tension are neglected. 

The only body force considered to be acting on the fluid is that due to gravity, 
and it acts along the negative z-axis as shown in Figure i. These assump- 
tions are justifiable [9] in the study of many fluid oscillation problems, such 
as the sloshing of the propellant in the tanks of a space vehicle during boost 
flight. 

The above assumptions ensure the existence of a velocity potential 
ip(r,6, z;t) for the fluid [10] . Thus, the velocity field can be expressed as 

T = V ip in V , (2.1) 

and the continuity equation, in terms of the velocity potential, becomes 
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Laplace’s equation 


V 2 ip = 0 in V (2.2) 

Thus, the equation of motion of the fluid is Laplace's equation, a linear 
differential equation about which a great deal is known* 

The boundary conditions that the solutions of (2.2) must satisfy are 
based on the assumptions that any particle on the boundary of the fluid remains 
on the boundary and that the pressure on the free surface is known. To 
determine mathematical Expressions for the boundary conditions, let the 
entire surface of the fluid, both fixed and free, be expressed as 

6(r,0,z;t) = 0 (2.3) 

Taking the material derivative of (2. 3) yields 

V6 • Vi/> + S = 0 . (2.4) 

b 

Noting that V <5 • is the velocity normal to the surface S , 
equation (2.4) becomes 

dip/dn + <5 = 0 on S (2.5) 

Since <5=0 on the fluid surface in contact with the container, the 

b 

boundary condition on F becomes 

dip/dn = 0 on F (2.6) 
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Letting ri(r,6;t) be the wave height of the free surface, 8 
may be expressed as 

8 = z - T) on S (2.7) 

Then substituting (2. 7) in (2. 5) , the boundary condition on the free 
surface becomes 

^r + ^e 77 / 1,2 ~ + v t = 0 ° n s • ( 2 -8) 

It is necessary to obtain another boundary condition for the free 
surface since the unknown wave height rj has been introduced into the 
problem. The additional boundary condition comes from the previously 
mentioned assumption that the pressure on the free surface is known. A 
mathematical expression for this boundary comes from Bernouilli's law 

+ (V^ • V ^)/2 + (P - P 0 )/p + gz = 0 , (2.9) 

where P 0 is the pressure at the free surface . On the free surface S , 
equation (2.9) becomes 

+ (V# • V^)/2 + rjg = 0 on S (2.10) 

The differential equation and appropriate boundary conditions that 
govern the motion of a fluid undergoing free oscillations in a rigid, partially 
filled container have been developed. These governing equations are 
summarized as follows: 
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HS- 

t! 

O 

in V 

(2.11a) 

dip/dn = 0 

on F 

(2.11b) 

V>r + W r2 - K + ”t ' 0 

on S 

(2.11c) 

ip + ( V ip • V ip)/ 2 + rjg = 0 

on S 

(2. lid) 


Equations (2. li) represent a boundary value problem in which the 
differential equation is linear, and the boundary conditions consist of one 
linear boundary condition and two nonlinear boundary conditions. Not only 
are equations (2. lie) and (2. lid) nonlinear, but they must be satisfied on a 
moving surface. Because of the difficulty of solving such a boundary value 
problem, very little progress has been made without making simplifying 
assumptions of one sort or another. 

The assumption is made in this analysis that the wave height r\ and 
the velocity potential tp can be expanded in a power series with respect to a 
small positive parameter e , which is proportional to the wave height of 
linear solution and becomes evident later in the analysis. 


ip = &P 1 + e 2 ip 2 + e 3 ip 3 + . . . (2.12) 

T) = et] 1 + e 2 T) 2 + e 3 rj 3 + (2.13) 


Equations (2.11a) and (2.11b), through the use of (2.12), become 
V 2 ^ =0 in V (2.14) 

dtp^/dn =0 on F (2.15) 
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Since boundary conditions (2. lie) and (2. lid) require evaluating the 
velocity potential and its derivatives on the free surface (z = rj) , it is 
assumed that these functions can be expanded in a Taylor's series about the 
undisturbed free surface (that is, z = 0) . Thus 

$(r, 0, 7j;t) = ip( r,0,O;t) + ip ( r,0,O;t)r } 

z 

+ $ (r, 0, 0;t)77 2 /2 + (2.16) 

zz 

Substituting equations (2.12) and (2.13) in (2.16), yields 

</>(r, 0, ? 7 ;t) = ef^r, 0, 0; t) ] + e 2 [^ 2 (r, 0, 0;t) + ip 1 (r, 0, 0; t) r? 1 ] 

z 

+ e 3 [^ 8 (r,0, 0; t) + ^ 2 (r, 0, 0; t) rj 1 + ^(r, 0, 0;t)rj 2 

z z 

+ (r, 0, 0;t) T] i ri i ] + 0(e 4 ) (2.17) 


Since the derivatives of ip have similar expansions, equation (2.11c) 
becomes 


e(7? t “ + e2(r, t " K _ O* + + 


1 .,1 


,l~i 


, 1 /_ 2 ■ 


+ e3(T? t “ “ O' - Kz^ - Kzz^ + KK + n rK 

+ V^V 1 + i? 0 2 ^Vr 2 + Tj^^/r 2 + ^^Vr 2 ) 


+ 0(e 4 ) = 0 


on S 0 


(2.18) 
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and equation (2. ild) becomes 


+ i^g) + e 2 {^ 2 + + gv 2 

+ [ ( 4>p 2 + ( $*) 2 /r 2 + ( ipp 2 ] / 2 } 

+ e 3 (^ 3 + ip 1 ^ 2 + ip 2 r } 1 + ^ 1 r] 1 r ] 1 + grj 3 + ip i ip z + 77* 

1/ tz iz zzt i* it i* rz 

+ ^i/r 2 + tyty 1 r? 1 / r 2 + i/) 1 ?/) 2 + tpi-ip 1 r) 1 ) 

t/C/ u u Z ZZ Z ZZ 


+ 0(e 4 ) = 0 

on S 0 

(2.19) 

Equating the coefficients of e 11 equal to zero, 

the equations to be 

solved for the first-order become 



vV = o 

in V 

(2.20a) 

9 !p i /dn = 0 

on F 

(2.20b) 

V, 1 ~ ip 1 = o 

t Z 

on S 0 

(2.20c) 

ipl + grj 1 = 0 

on S 0 

( 2. 20d) 

For the second-order, the equations become 



V 2 ip 2 = 0 

in V 

(2.21a) 

dip 2 /dn = 0 

on F 

(2.21b) 

^ = KvP 1 - 

+ g T? 2 = - - [( ^) 2 + ( 2 /r 2 

on S 0 

(2.21c) 

+ ( ^) 2 ] 

on S 0 

(2. 21d) 
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And, for the third-order, the equations become 


V 2 ^ 3 = 0 
8^ 3 /8n = 0 


in V (2.22a) 

on F (2.22b) 


T7 3 - jp 2 = ip 1 T7 2 + ip 2 ri l + ip* ri i r] i - ri 2 ip 1 - v 1 ^ 2 
t t zz zz zzz r r r r 

- - Vgtg/r 2 - V**l>*/r* ~ Vq $q z V 1 /* 2 

on Sj 


(2.22c) 


^ t 3 + gv 3 


- C”’ - K/ ' *zzt nW 


' -V* • -v» ' -v» ' -v» r» ' 


*X /r2 ’ ^9z’ ,1/rJ ' ‘ 


on S 0 


(2. 22d) 


Equations (2.20) through (2.22), in principle, furnish a means of 
calculating successively the coefficients of the series (2. 12) and (2. 13) . 
Equations (2.20) lead to solutions for ip 1 and r j 1 . Once i ^ 1 and t] 1 are 
determined, they can be inserted in equations ( 2 . 21 ) to permit solutions for 
ip 2 and T ] 2 • Thus, this procedure can be applied until as many terms in the 
series for ip and 17 as desired have been calculated. 

Even though the procedure described above will yield a solution, there 
remain two boundary conditions on the free surface and two unknowns in the 
problem to determine. Therefore, it is a computational advantage to combine 
the two free surface boundary conditions so that 77 is eliminated, leaving a 
boundary value problem with one unknown to determine. This is simple to do 
once equations (2.20d), (2.21d), and (2.22d) have been differentiated with 
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respect to time. It must be remembered that combining the two free surface 
boundary equations in this manner eliminates any constant term which may 
exist in tj 1 , r\ 2 , and r\ z . These constant terms will be dealt with later. 


The boundary value problem to be solved, after the algebriac manipulations 
have been performed, is for the first-order, 



vV= o 

in V 

(2.23a) 


0^/0 n = 0 

on F 

(2.23b) 


< + < = 0 

on S 0 

(2.23c) 

and 

II 

1 

on S 0 

(2.23d) 

For the second-order, the problem is 




vV 2 = o 

in V 

(2. 24a) 


dip 2 / 9n = 0 

on F 

(2.24b) 


♦A + ‘ + ' 2 *r*rt ‘ 

2 W/ r * 


and 

- 2 KKt 

on S 0 

(2.24c) 


r? = - # t Vg - ♦ t ^ t 1 /g ! - [ ( # r ‘) 2 + ( ♦ *) 2 /r 

2 + (^) 2 ]/(2g) 
z 



on S 0 

( 2 . 24d) 


And, for the third-order, the problem is 
V 2 ip 3 = 0 


in V (2.25a) 
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dip s /dn = 0 


on F 


(2.25b) 


+ < 


2 *r*rt - - 2^/r* ' ^eV 1 * 


- 2 *z # zt ' 2# z Kt + /g + 

+ *zzH 1 + Kfl ' ( >2 ^rr " ‘♦(W/ 1 * 

+ < 5 e ‘ t^ji/r 3 - y*y z - wy*y * 2 

- 2 ^sW r! + ^ttz^r^* 2 ** + ( ^ * ! ''ttz /,< 2gr2) 

+ f t J z < j 1 ^ )V ( 2g) + ? z l z < *! ) ! /2 + ^ z (# 0 l )V(2r 2 ) 

+ ^ z ( ^)V2 + w t xy/g + y*y z y 

+ 2 WA‘t /(gr2) + 2 *tWz/ (gr,) + 2 *t*L' l 'zt /g 
* 2 *X*Lt /g - ■ ^Vt 1 ^/ 8 


- (* t ‘>^ zz /(2 g ) 


i \2,/, i 


on S 0 


(2.25c) 


and 


7? 3 = ~ fyVg * ^^ t 2 / g 2 - ( ^ z ) 2 ^/g 3 


+ ^ t 2[(^ )2 + (^) 2 / r2 + (^ )2] / (2g2) + ^ tZ ^ t ^ 2 

+ ^ z 1 zt ( ^ )Vg 2 + ^ 2 + ^^rz^t /g2 0n s ° 


( 2. 25d) 


The constant terms that were removed from rj 1 , rf , and rj 2 in the 
process of combining the free surface boundary equations may be lumped 
together in one constant, 77 0 . The wave height equation then becomes 
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V 


( 2 . 26 ) 


= 7)° + erj 1 + e 2 r? 2 + e 3 i 7 3 + . . . on S 0 

The constant rj° is evaluated from the condition that the volume of 
fluid remains constant and can be expressed mathematically as 

Jr ?dS 0 = 0 . (2.27) 

So 

The equations governing the nonlinear oscillations of a fluid in a 
partially filled container have been developed. 
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LINEAR SOLUTION 


It is common in many fluid oscillation problems to retain only the 
first-order terms in the governing equations . The problem is then known as 
a linear fluid oscillation problem. This linearization is accomplished by 
truncating the series expansions of ip and rj after the first term; therefore, 
the governing equations are the first-order equations expressed by (2. 23) . 
These equations are restated as 


vV = o 

in V 

(3.1a) 

dip 1 /da = 0 

on F 

(3.1b) 

+ S K = 0 

on S 0 

(3.1c) 


and 

V 1 = on S 0 (3. id) 

Smce this analysis is concerned only with fluid motion that is 
harmonic with time, the time dependence can be removed from the equations 
by the transformation 

^ = 0 (r,0, z)e lwt , (3.2) 
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where co is the frequency of oscillation. Then equations ( 3. 1) become 


V 2 0 =0 in V (3.3a) 

d<p/dn = 0 on F (3.3b) 

<p = oj 2 cp/g on S 0 . (3.3c) 

z 

The boundary value problem (3.3) can also be expressed as an 
extremum problem. According to Hamilton's principle 

61 = 0 , (3.4) 

where 61 is an isochronous variation of 

I = f( KE - PE)dt . (3.5) 

t 

The kinetic and potential energy are expressed as 

KE = (p/2) f7 • 7dV (3.6) 

V 

PE = (pg/2) f (V 1 ) 2 dS 0 , (3.7) 

So 

where p is the fluid density. 

If the kinetic and potential energy are written in terms of the velocity 
potential, equation (3.5) becomes 

I = /[p/2/(V^ • V ip 1 ) dV - (p/2g) / ( # 1 ) 2 dS 0 ] dt . (3.8) 

tv S 0 
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By integrating over t from 0 to 2ir/co and omitting a nonessential 


factor, I can be expressed as 

I = /(V0 • V <p) dV - (cu 2 /g) / ( 0) 2 dS 0 . (3.9) 

V S 0 

Lawrence, Wang, and Reddy [11] have shown that the solutions of 
equations (3.3) are also the solutions of the extremum problem for the 
integral I as expressed by equation (3.9) . 

The approach taken in this analysis is to determine an approximate 
solution to the extremum problem rather than to solve the boundary value 
equations. Trefftz’s method [12] will be used to construct this approximate 
solution. This method begins by choosing a complete sequence of linearly 
independent functions, u.. (r, 0,z), each of which satisfies the differential 

IK 

equation. The approximate solution is expressed as a series of the 

coordinate functions u.. 

lk 

^ = Hvik ■ (3 - 10) 

i k 

where the coefficients a., are determined from the condition 

lk 

1(0 - $) = minimum (3.11) 


A set of necessary conditions are 
81(0 - $)/da = 0 


(3. 12) 
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Through the use of equations (3. 9) and (3. 10) , equation (3. 12) 
becomes 




Vu ik ) dV “ (C ° 2 g) 


( u jiV ds »J 


a = 0 . (3.13) 


The volume integral in (3. 13) can be changed to a surface integral 
using Green's Theorem. Applying this theorem to the volume integral and 
noting that u^ and u.^ satisfy Laplace's equation and that the dot product 
obeys the commutative law, the volume integral becomes 

f (Vu.. • Vu.,)dV = f ( u.. 9u.. /3n + u.. 3u,./9n)dF 
J jl ik 9. J il ik lk il 


V 


2 Jl ^ 


+ 2 / (u jl 8u u/ 8n + u ik au i/ 9n)dS ° 

So 


ik jl 

... 3u , 
ik jl 


(3.14) 


Equation (3. 13) can now be written as 


i 1 1 F ' 


... 3u.. /8n 
jl ik 


u.. 8u. T /8n) dF 
ik jl 


+ 



8u.. /8n + 
ik 


U ik 9u j^ n ) dS 0 



(3.15) 


Since this analysis is concerned with axisymmetric containers, it is 
advantageous to consider the surface F to be generated by rotating a curve 
y about the axis of symmetry. 
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Let y be expressed in parametric form 

r = r(o!) 
z = z(o;) 

where 

aj < ce ^ a 2 

and 

r(o! 1 ) = 0 z( ex. i) = -H 

r(a 2 ) = R z(a 2 ) = 0 

where R is the radius of the cbntainer at the free surface, and H is the 
height of the container at the axis of symmetry. 

The derivatives normal to F and differential element of area on F 
then become 

8( )/8n = [z’8( )/8r - r'8( )/8z]/V(r') 2 + (z’) 2 (3.16) 

dF = W (r*) 2 + (z 1 ) 2 dedo: (3.17) 

where 

r 1 = dr(o!)/dQ! z 1 = dz(a)/do! 

Using equations (3. 16) and (3.17) and realizing that the surface S 0 
is a circle, we can express equation (3. 15) as 
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{ 2ir (*2 

f I [u (z’au. k /ar - r'au /az) 

o J 

+ u.. (z' 9u.,/9r - r'au.Vaz)] rdcude 
ik jl jl 

27 r R 

+ f f [u.. 9u.,/9z + u„ 9u../9z ] rdrd0 
J 0 J 0 Jl ik ik jl 

2tt R | 

- (2c^/g) J / (u^u^) rdrdej = 0 (3.18) 

The chief difficulty encountered in solving equations (3. 18) is 
selecting the coordinate functions since there are no general recommendations 
available . It is known from the general theory that if the system of coordinate 
functions is complete and it satisfies the differential equation, the minimizing 
sequence constructed by Trefftz's method will converge to the exact solution 
[12] . It is therefore advisable to select u ik as eigenfunctions of some volume 
that contains the given volume but has a simpler shape. Therefore, in this 
analysis, the coordinate functions that will be used are the eigenfunctions of 
the fluid in a cylindrical container whose cross section equals the largest 
cross section of the given container. Thus 

U ik = J k (i8 ki r) cosh[ / 3 k i (Z + H c ) ] cos(k0)/cosh(/3 k .H c ) (3.19) 

where B are the roots of dJ (/3 r)/dr I = 0 ; R is the maximum radius 

kl k r=R C 

c 

of the axisymmetric container; and is the depth of the fluid in the 
axisymmetric container. 
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Equation (3.19) can be rewritten as 


U ik = ^ik cos(k0) * (3.20) 

where 

V = J k ( ^ki r)0O8h ^kl (z + H o ,1/cOSh( ^ki H c ) 

Inserting the coordinate functions expressed by (3.20) into (3.18) 
and integrating with respect to 9 , yields 


OO OO 

z z 


j=l k=0 


+- 


f [u., (z'8u., /dr - r'8u../8z) 

J lk lk lk 

O'! 

+ u (z 1 9u.,/9r - r'8u /8z) ] rdcn 

IK jK JK 

f < V V 8z + s ik a V 3z)rdr 

R \ 


- (2coVg)/(u. k u. k )rdrJa jk = 0 


(3.21) 


where a nonessential factor has been dropped and the orthogonal 'properties 
of the cosine functions have been used. 

Equation (3.21) represents an eigenvalue problem for each k . 
Since equation (3.21) is symmetric with respect to i and j , the eigen- 
values, 2a> 2 /g, are real. Solving equation (3. 21) by one of several 
standard means for each value of k , the eigenvalues and vectors are 
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2\l./g ~ ei S envalue 


a — jth element of the eigenvector 

jkl 


Thus the solutions of (3.3) are 


0 k l = .2 a jk l J k O k j r) COSh[/3 ki ( z + H n )] cos(k0) /cosh( ^ ki H c> 


kj' 


kj c' 

(3.22) 


and the frequency of oscillation is 


W kl = ^kl 


(3.23) 


Equations (3.22) and (3.23) are the eigenfunctions and the frequency 
of oscillation of a fluid in an axisymmetric container. 

The eigenfunctions have the following important orthogonality 
relationship which will be needed later: 


/ V s * = \i s ii < 3 - 2 ' 

b o 

This relationship is easily obtained by using the orthogonal 
properties of the cosine functions and the eigenvectors. 

A means has been presented by which the linear fluid oscillation 
problem can be solved for the case of a fluid contained in an axisymmetric 
vessel of general shape. 
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NONLINEAR SOLUTION 


The equations governing the nonlinear fluid oscillation problem 
were developed previously. The velocity potential and wave height were i 
assumed to be expandable in a power series of a small positive parameter, 
and the governing equations were developed for the terms in the expansion 
for the velocity potential through the third -order term. The method by which 
the nonlinear equations will be solved is that of Krylov and Bogoliubov [ 8] . 
This method assumes that velocity potential can be expressed as 


,.(i> 


OO 

V\ 


(i) 


= y. Y. a ; <p 

p=0 m=i pm pm 


( 4 . 1 ) 


It is obvious that the assumed form of satisfies the differential 


equation and the boundary condition on F , since are solutions of the 

linear problem. Thus, the problem at hand is to evaluate so that the 

free surface boundary condition is satisfied. Since we are seeking periodic 
solutions, is periodic with time; therefore, it is advantageous to effect 

a change of variable so that is periodic with a period of 2n . Thus let 
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T 


(4.2) 


t 

= /(co + eB 1 + e 2 B 2 + . . . ) dt 
0 

where co is the linear frequency. 

Therefore, the instantaneous frequency f is obtained by differentiating 
equation (4.2) 

r = co + eB 1 + e 2 B 2 + . . . (4.3) 

In general, the terms B^ are functions of the wave amplitude. 

Thus, for the sake of completeness, let equation (4.3) be written as 

r = co + eB^K) + e 2 B 2 (K) + . . . (4.4) 


where K is an amplitude parameter, and its rate of change is a function of 
the amplitude . Thus K can be defined by 


K = eD 1 ( K) + e 2 D 2 (K) + . . . 


(4.5) 


Equation (4. 1) can now be written as 


* 


(i) 


oo oo 


Z 


p=0 m=l 


pm 


pm 


where the parameters r and K are defined by 


(4.6) 


r = co + cB^K) + e 2 B 2 (K) + . . . (4.7) 

K = eD 1 (K) + e 2 D 2 (K) + (4.8) 
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To uniquely define the coefficients in equations (4. 7) and (4. 8) it is 
necessary to impose additional restrictions on them. Mathematically, these 
conditions can be expressed as 

2tt 

J J ipK cost d S 0 dr = eK 2 7r (4.9) 

0 S 0 

2tt 

f f ipK sinr <pu dS 0 dT = 0 (4.10) 

0 S 0 

From a physical point of view, the imposition of these conditions is equivalent 
to selecting K as the full amplitude of the first harmonic of the oscillation. 

Noting that ip = ip[r,6 , z,r(t) ,K(t) ] , the time derivatives of the 
velocity potential become 

\ = K' r + (4,11) 

*tt = ( Vt* + • (4 ’ 12) 

Substituting the expansions for f , K, and ip in equations (4. 11) 
and (4.12) yields 

= eiioip 1 ) + e 2 (o>^ 2 + B 1 ^ 1 + Dto 

t T T T K 

+ £ 3 (coi/> 3 + Bty 2 + B v + dVJ + D 2 ^) + 0(e 4 ) 

T T T K & 

(4.13) 
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and 


!/» = e(c o 2 !/> ) + e 2 (u ) Z 4> 2 + 2wBV + 2uD^J: ) 

tt TT TT TT IVT 

+ e 3 [co 2 ^ 3 + D j B t ^ P 1 + (B 1 ) 2 ip i + 2coB z ip 1 

TT K T TT TT 

+ 2wB^ 2 + 2coD^ 2 + 2B 1 DVJ + 2coDVj 
TT Kt Kt Kt 

+ (D 1 ) 2 ^ + + 0(e4) (4 - 14) 


Substituting the expansions for ^ and ip in equations (4. 13) and 
(4.14) and equating coefficients of e 1 , yields 


ip} = ojip 1 

t T 


*; = 


r = 


co tp Z + B 1 ^ 1 + D 1 !/)* 
T T K 


0) ip 3 + bV + B^ 1 + + D z ii; 

T T T K K 


52,/, i 


\Uh 2 


I 2 ,/, i 


(4.15) 

(4.16) 

(4.17) 


*£ = “K 


(4.18) 


ip z = co 2 !/' 2 + 2coB 1 i/> 1 + 2co D 1 !/) * 

tt TT TT Kt 


(4.19) 


i /' 3 = CO 2 !/' 3 + D*B* !/>* + 2a>BV + ( B 1 ) 2 J/' 1 + 2coBV 

tt TT K T TT TT TT 

+ 2wD 1 p z + 2B 1 D 1 i/'J + 2coD 2 i/'l + (D 1 ) 2 ^ 

I\T IVT IVT rvJv 


+ d ‘ d k^ 


(4.20) 


Through the use of equations (4. 15) through (4. 20) , the boundary 
condition equations (2.23c), (2.24c), and (2.25c) become, for the first-order 


27 



«\\ + = ° 


and, for the second -order 


“% T + g^ z 


2B 1 c dp 1 - 2coD 1 i/>j - 2c oip 1 ip 1 
tt Kt t r r 


-2toip 1 „^„Vr 2 - 2co ipj ip 1 + (*>H 1 ^Vg 
T0 0 Tz Z TTZ T 


+ CO 1 
T ZZ 


(4.22) 


and, finally, for the third -order 


to 2 ip 3 + gip s 

TT Z 


- B^D 1 ^ 1 - 2coBV - (B 1 ) 2 ^ 1 - 2toB i ip 2 

XV T TT TT TT 


-2o.D‘^ t - 2b’D 1 ^ t - 2wD 2 ^ t - (d') ! ^ k 

■ D ‘ D K # K ' 2B ‘*i* r ‘ ' 2 "* T Vr‘ " 


2w^ 1 i p 2 - 2B i ip i ip^/x 2 - 2toip 2 4>}/r 2 

>r r r V0 V t0 6 ' 


2D 1 ip*ip}/T z - 2toip 1 n ip 2 /x' i - 2B i ip i ip 1 
K0 r 0 r T0 r 0 tz”z 


-2 GJip 2 Ip 1 - 2D7 T i tp 1 - 2coip 1 jp 2 

TZ Z Kz Z TZ Z 


+ 2o) 2 B 1 !/> 1 ^Vg + Cl) 3 !/' 2 ^ Vg + 2ci) 2 D 1 !/' T j !/>Vg 

TTZ T TTZ T KTZ T 


+ Ci ) 3 !/ 1 1 ^> 2 /g + Ci) 2 B 1 !/» 1 !/)Vg + Ci) 2 D 1 !/ 1 ip^/g 

TTZ T TTZ T TTZ XV 


+ Ci jip 1 ^ 2 + to ip 2 ^ 1 + B 1 !/ 1 ^ 1 + D 1 ^*^ 1 

T ZZ T ZZ T ZZ K ZZ 


(iph 2 ^ 1 - (i p 1 ) 2 ip i /r 4 - (ip^) 2 ip^ + (ip i ) 2 ip 1 /x 3 

r’ rr V V W V V zz 'V V 


2 ^r‘z ■ - KWJ 1 


1,1, 1,1, 1 /JL _ ,1, 1,1, 1,1, 1 /-v.2 


+ co 2 tp 1 (tp* )7( 2g) + o) 2 ( i/ 1 ) V /( 2gr 2 ) 

TTZ r 0 TTZ 
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+ coVj0 2 /(2g) + (^)Vi-/ 2 + <0 2i tV< 2r2 ) 

TTZ Z X ZZ (7 ZZ 


+ (ip 1 ) 2 !/) 1 /2 + 2co i !p 1 ip 1 ip 1 / g + 2co 2 ip 1 ip i ip 1 /g 
z zz r rz rr r rrzr 


+ 2oj2 ^6zKo^ gr2) + 2co2 ^T^ 1 ^T0z /(gr2) 


+ 2co 2 ip 1 tp 1 ip 1 /g + 2<jj 2 ip 1 ip 1 ip 1 /g 
T ZZ TZ T Z TZZ 


- OJ 4 ^ ^ 1 ^ 1 /g 2 - aj 2 ip 1 i/) 1 !/) 1 /g 

T TZ TTZ ° TZ T ZZ 

- co 4 ( 2 !/' r 1 Tzz /( 2g 2 ) - co 2 ( V z 1 zz /( 2g) 

(4.23) 

As expressed by (4.6), ^ ^ satisfies the differential equation and 
the boundary condition at the container wall F . Thus:, if the terms 
are chosen so that the boundary condition oh the free slurface equations (4. 21) , 
( 4. 22) , and ( 4. 23) are satisfied, then the nonlinear fluid oscillation problem 
has been solved through the third -order terms. 

Since this analysis is to determine a solution near the first non- 
axisymmetric mode of oscillation, c o is taken to be defined by 


CO = COil 


(4.24) 


First-Order Solution 

It is clear from the linear solution shown previously that equation 
(4.21) is satisfied by choosing 
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ip 1 = K cos (t) <pu 


(4.25) 


The wave height rj 1 can be obtained from equation (2.23d) after equations 
(4.15) and (4.25) have been used. 


Second-Order Solution 

To determine an expression for A 2 so that the second-order free 

pm 

surface boundary condition is satisfied, equations (4.6) and (4.24) are 
substituted in equation ( 4. 22) to yield 


oo OO 


2 2 


{[w 2 (A 2 ) 

pm tt 


+ 


(JO 2 A 2 ] $ 

pm pm pm 


cos(p0) } 


= sin(2r){a;K 2 [4 11 ) r ] 2 /2 + coK 2 ( $ n ) 2 /( 2r 2 ) + a> 5 K 2 ($ n ) 2 /( 2g 2 ) } 

+ cos (0)[2B 1 O)K^ 11 cos(t) + 2a>D 1 K$ 11 sin(r) ] 

+ cos (20) sin(2T){wK 2 [($ 11 ) r ] 2 /2 + co 2 K($ 11 ) 2 /(2r 2 ) 


+ u) 5 K 2 ($ 11 ) 2 /(2g 2 )} , (4.26) 


where the identities 


cos 2 (p0) 

= [1 + cos(2p0)]/2 

(4.27) 

sin 2 (p0) 

= [1 - cos(2p0)]/2 

(4.28) 

sin(2T) 

= 2 sin(T) cos(t) 

(4.29) 
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and the expressions 


cb = d> cos(p0) 
pm pm r ' 


(<p ) = c o 2 <p /g 

pm z pm ^pm 


(4.30) 


(4.31) 


have been used. 

The two summations in equation ( 4. 26) can be removed by using the 

orthogonality relationship (3. 24) after the equation has been multiplied by 

<p and integrated over the surface S 0 . The result is 
qn ? 

w 2 (A 2 ) + A 2 = 2B 1 oo K cos(r) <5 , 6 + 2ojD 1 sin(r)6 , 6 . 

qn tt qn qn qi nl ' ' ql nl 

+ tto;K 2 sin(2r)(I 1 + I 2 + co 4 I 3 /g 2 )6 

qn qn qn qO 

+ 7ro;K 2 sin(2T) (I 1 - l 2 

qn qn 


ca 4 I 3 /g 2 )6 /2 

qn & q2 


( 4. 32) 


where 


I 1 

qn 



I 2 

qn 



I 3 

qn 


R 



($n) 2 rdr 


(4.33) 


(4.34) 


(4.35) 
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Substituting the Fourier expansion of 


A 2 

qn 


|a° + 

2 qn 



+ 



sin(i t) ] 


in equation ( 4. 32) , it becomes 


(4. 36) 


ct> 2 a°/2 + /’"'[(c 2 -f 2 c 2 ) a^ cos(ir) + (c 2 -i 2 c 2 )b^ sin(i r) ] 

qn qn 1 qn qn qn qn 

= 6 6 2c [KB 1 cos(t) + D 1 sin(T)] + ttcK 2 6 .(I 1 + I 2 

qi ni qO qn qn 

+ o; i I ( ^ n /g 2 )sin(2T) 

(i/2)7rcoK 2 6 q2 (I^ - I* n + w 4 I^ n /g 2 )sin(2T) (4.37) 


The Fourier coefficients of the expansion for A 2 can be determined 

qn 

from equation ( 4. 37) by equating the coefficients of cos(fr) and sin(f r) 
after the conditions expressed by equations (4. 9) and (4. 10) have been 
imposed and the assumption made that w 2 ' n # H 2 u ; 2 . This assumption will be 
examined later. Thus, 


B 1 

D 1 


0 

0 


™ K2(I 0n + I 0n + AV^Kn ' 4 " 2 > 

(i/2) ir oj 2 K 2 (I 2n - Ij n + ^ 4 I 2 3 n /g 2 )/(co 2 2 n - 4c 2 ) 


(4.38) 

(4.39) 

(4.40) 

(4.41) 
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0 


(4,42) 


i* 2 (4.43) 

q # 0 or q + 2 (4. 44) 

When equations (4.38) through (4.44) are used, the second-order term 
in the expansion of the velocity potential becomes 

OO 

^ * sto < 2T > 2 (b to*0n + b L+2n> (4 ' 45) 

n=i 

The wave height rf can be evaluated from equation (2.24d) after equations 
(4.15), (4.16), (4.25), and (4.45) have been used. 

Third-Order Solution 

The third-order solution is developed in a manner similar to that 
used for the second -order solution. Equations (4.6), (4.25), and (4.45) 
are substituted in the third-order free surface boundary condition (4.23) to 
yield 


SL 

qn 

i 

'qn 


a 


b 


= 0 


.2 _ 


qn 


= 0 


V y {[oj 2 (A 3 ) + cu 2 A 3 ] $ cos(p0)} 

U U 1 1 ' pm ' tt pm pm J ^pm 


p=0 m=i 


= cos(t) cos(0){2coKB 2 ^ 1 - wK(^ 11 ) r b oj ( ^ 0 j ) r 


- ( 1/2) a) K( (j>u) r E^b|j ( v* 2 j) r 


33 



oo 


- (l/4)a ) K^ 1 ^b 2 2 . u *.^ 2 .|/g ! 

- (9/16)K 3 [(0ii) r ] 2 (0ii) rr + O/iejKVt^Vg 4 

- (3/8)KV[(^ 11 ) r ] 2 0 11 /g 2 + (3/16)K 3 ($ 11 ) 3 /r 4 

+ (3/i6)K 3 ($ 11 ) 2 ($ 11 ) r /r 3 - ( 3/ 8) K 3 [( $n) p ] 2 $u/ r 2 

- (i/8)K 3 o; 4 ($ 11 ) 3 /(g 2 r 2 ) ■ 

j 

+ cos(r) cos(30)| - (i/2)coK($ 11 ) r _E b 2 j ( ^ 2 j ) r 

+ (i/2) U ^^|b| jW * j $ 2 ^ 

- ( 1/4) co 4i^J> 2 .co 4 . $^j/g 2 

-(3/l6)K 3 [(J 11 ) j J 2 (J 11 ) rr 
+ (1/16) K 3 co 8 (^n) 3 /g 4 

- (i/8)K 3 u 4 [($ ii ) r ] 2 $ n /g 2 - (3/16)K 8 ($ u ) 8 /r 4 
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- (3/i6)K 8 ($ 11 ) 2 ($ 11 ) r /r 8 
+ (3/8)K 3 [($ 11 ) r ] 2 $n/r 2 

+ (i/8)K 3 w 4 ($ 11 ) 3 /(g 2 r 2 ) J 

+ cos(3t) cos( 6 ) j-3coK($ n ) r J) b Qj ( ^oj ) r 

OO 

-(S/2)($„) r Zb°.($ 2 .) r 

- I J^/r 2 

- (5/2) 

- 5a,*K$ 11 ^|b^a, 2 .$ 0} ^ 


*( 1 / 2 )^$^^.$^ 

+ ( 1/4) a.4^2 b 2 .^ 1 . 

- (3/i6)K 3 [($ 11 ) r ] 2 (J 11 ) rr 

- (15/l6)K 3 co 8 (J 11 ) 3 /g 4 

- (9/8)K 3 co 4 J 11 [(J 11 ) r ] 2 /g 2 + (i/i6)K 3 (0 11 ) 3 /r 4 
+ (i/16)K 3 (J 11 ) 2 (J 11 ) r /r 3 
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- (i/8)K 3 [(J 11 ) r ] 2 Jn/r 2 

- ( 3/8) K s w 4 ( $„)*/( r 2 g*) 


+ cos(3r) cos(30) | - (3/2)cuK<J 11 ) r 2 b 2j ^ 2 ]^ 


ScoK^/f b * $\/t 


- (5/2) W 3 K0 1 |^b 2 .co 2 2 .0 2 y g 2 


-(l/l6)K 3 [(^ 11 ) r ] 2 (0 11 ) rr 

- (5/16)K 3 w 8 (^ 11 ) 3 /g 4 

- (3/8)u i K 3 [($ n ) r ?$ n /g 2 

- (l/l^K^^Vr 4 - (i/i6)K 3 (0 11 ) 2 (^ 11 ) r /r 3 
-i- (i/8)K 3 [(^ 11 ) r ] 2 0 11 /r 2 

(3/8)KV({ tl ) 3 /(g 2 r 2 ) | 

+ sin(r) cos(0) (2coD 2 ^ 11 ) . (4.46) 


The two summations in equation (4.46) can be removed by using the 
orthogonality relationship (3.24) after the equation has been multiplied by 
and integrated over the surface S 0 . The result is 
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coKI 6 

qn 


co 2 (A 3 ) + co 2 A 3 

qn tt qn qn 

=, cos(t){2coKB 2 6 ,5 , + 8 ,tt[-coKI 4 - (i/2)coKI 5 - 
qi ni qi qn qn 

+ (i/2)c0 3 KI 7 /g 2 + b 3 KlVg 2 


- (i/2)coKL^/g 2 - ( 1/4) co K 1^/ g 2 

-(e/iGJK 3 ! 11 + (3/16) K 3 co 8 1 13 /g 4 
qn qn 

- (3/8)K 3 co 4 I 15 /g 2 + (3/16)K 3 I 12 

qn ° qn 

+ (3/16)K 3 I 14 - (3/8)K 3 I le 
qn qn 

- (1/8) K 3 co 4 1^/ g 2 ] 

+ 6 0 7r [ - ( 1/2) co KI 5 + coKI 6 + (i/2)co 3 KI 7 /g 2 

q3 7 qn qn qn ° 

- (i/4)coKI^/g 2 - (3/16)K 3 I^ 

+ (l/16)K 3 co 8 I 13 /g 4 - (l/8)K 3 co 4 I 15 /g 2 
qn qn 

- ( 3/16) K 3 1 12 - ( 3/16) K 3 1 14 + ( 3/8) K 3 1 16 

qn qn qn 

+ (l/8)K 3 co 4 I^ 7 /g 2 ]}. 

+ cos( 3r) { 8 7r [ - 3co K 1 4 - (3/2)coKI 5 - 3coKI 6 

qi qn qn qn 

- (5/2)co 3 KlJ n /g 2 - 5co 3 KlVg 2 

+ ( 1/2) co K I^/g 2 + ( 1/ 4) co K I^/g 2 

-(S/iejK 3 ! 11 - ( 15/16) K 3 co 8 1 13 /g 4 
qn qn 5 

- (9/8)K 3 co 4 I^/g 2 + (i/16)K 3 lJ 2 
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where 


+ ( 1/1 6) K 3 1 14 - ( 1/8) K 3 1 16 - ( 3/8) K 3 co 4 1 17 /g 2 ] 
qn qn qn & J 

+ 6 „7T[- (3/2)u)KI 5 + 3u) KI 6 - ( 5/2) co 3 KI 7 /g 2 

q3 qn qn qn 

+ (i/4)w KI^/g 2 - (i/16)K 3 I^ 

- (5/16)K 3 o) 8 I 13 /g 4 - ( 3/ 8 ) o) 4 K 3 1 15 /g 2 

qn & v qn 

- ( 1/1 6) K 3 1 12 - ( 1/1 6) K 3 1 14 + (1/8)K 3 I 16 

qn qn ' qn 

+ (3/8)KVl^/g 2 ]} 

+ sin(r) 2o)D 2 6 ql 6 nl , (4.47) 



(4.48) 


( 4. 49) 


(4. 50) 


(4.51) 


(4. 52) 
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(4.53) 


(4.54) 


(4. 55) 


(4. 56) 


(4.57) 


(4.58) 


(4.59) 


(4.60) 


(4.61) 
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Substituting the Fourier expansion of 

oo 

A 3 =. ( 1/2) c 0 + V [c^ cos(i t) + d 1 sin (It)] , (4.62) 

qn qn ^ 1 qn qn ' J ’ ' 

in equation (4.47) , the following expression is obtained: 

CO 

( 1/2) c 0 w 2 + y [(w 2 - iV) C i cos(fr) 

qn qn £-1 qn qn v 

+ (w 2 - i 2 co 2 )d^ sin(ir)] 

qn qn 

= cos(t){2wKB 2 6 ,S , + <5 ,7r[-coKI 4 - (i/2)coKI 5 - wKI 6 
ql nl ql qn qn qn 

+ (l/2)o) 3 KlJ n /g 2 + a; 3 KI 3 n /g 2 

- (l/2)coKI 9 /g 2 - ( 1/ 4) a; KI*V g 2 

qn qn 

- (e/lGJK 3 ! 11 + (3/ 16)K 3 co 8 I 13 /g 4 

qn qn & 

- (3/8)K 3 w 4 I 15 /g 2 + (3/16) K 3 I 12 

qn & qn 

+ (3/16) K 3 I 14 - (3/8) K 3 I le 
qn qn 

- ( 1/ 8) K 3 co 4 1 17 /g 2 ] 

7 qn & J 

+ 6 7r [ - (i/2)a;KI 5 + wKI 6 + ( 1/2) c^KI 7 /g 2 

q3 qn qn qn & 

- (l/4)co KI^/g 2 - (3/16)K 3 I^ 

+ (l/16)K 3 co 8 I 13 /g 4 - (l/SJKVlf/g 2 
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- (3/i6)K 3 I 12 - ( 3/ 1 6) K 3 1 14 + (3/8)K 3 I ie 

v qn qn qn 

+ (l/8)KVl"/g 2 ]} 

+ cos(3t){ 6 ql 7r[-3o;KI^ n - (3/2)coKI 3 n - 3coKI 3 n 

- (5/2)co 3 KI^/g 2 - 5co 3 KI 3 n /g 2 

+ (i/2) W KI 3 n /g 2 + (i/4)coKI^/g 2 

- (3/ 16) K 3 I q ^ " ( 15/16) K 3 u 8 l‘ 2 /g 4 

- (9/8)KVl^/g 2 + (1/16)K 3 I^ 

+ ( 1/ 16) K 8 1 14 - (i/8)K 3 I 16 - ( 3/ 8) K 3 w 4 I * V g 2 ] 

v qn qn qn 

+ 6 q3 7r[- (3/2)coKI^ + 3coKI« n - (5/2)Al^/g 2 

+ ( i/4) co KI^/g 2 - ( 1/16) K 8 l£ 

- (5/16) K 3 oo 8 I^/g 4 ' ( 3/8) co 4 K 3 I q Vg 2 

- ( 1/16) K 3 I q ^ - (l/16)K 3 l£ + (1/8)K 3 I^ 

+ ( 3/ 8) K 3 co 4 1 1 Vg 2 ] } 
qn 

+ sin(r) 2coD 2 6^6^ (4.63) 


The Fourier coefficients of the expansion for can be determined 

from equation ( 4. 63) by equating the coefficients of cos(4t) and sin(ir) 
after the conditions expressed by equations (4. 9) and (4. 10) have been 
imposed and the assumption made that co^ q * l 2 co 2 . This assumption will be 
examined later. Thus, 
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c J n = *[-«KI* n - (i/2)«KI | n - «KI« n + (i/WKI^/g 2 
+ w S K I« n /g* - ( i/2) co K I^/g 2 - ( i/4) co K l“/g* 
-O/iejK 8 !^ + (3/16)K 8 l£ + (S/iejKVl 18 /* 4 
+ (3/i6)K 8 I 1 “ - (3/8) KVl^/g 2 - (3/8)K 8 l“ 

- (l/SJKVl^/g^/Cco^ - co 2 ) 

n*i (4,64) 

°3n = *[-<i/2)coKlJ n + coKl| n + ( i/2) co 3 KlJ n /g 2 

- (i/4)coKI 3 ^/g 2 - (3/i6)K 3 I 3 “ - (3/16)K 8 l£ 

+ (i/l6)K 3 co 8 I 3 13 /g 4 - (3/16)K 8 I g “ - (i/8)K 3 co 4 I^/g 2 

+ (3/8)K 3 I 3 13 + (l/8)K 3 co 4 I 3 1 yg 2 ]/(co 3 2 n - co 2 ) 

° 1 (4.65) 

c 3 n = 7r[-3coKI 4 n - (3/2)a)KI 8 n - ScoKI^ - ( 5/2) co^I^/g 2 

- 5co 3 KI, 8 n /g 2 + ( i/2) co K I ®^/ g 2 + ( i/4) co K I 4 ® 

- (3/i6)K 8 I 4 ^ + (i/i6)K : 8 l£ - (lS/iejKVl^/g 4 

+ (l/i6)K 3 I, 14 - (9/8)K 3 co 4 I, 15 /g 2 - (i/8)K 3 I, 16 
in in In 

- (3/8)K 3 co 4 I 1 l y g 2 ]/(co 2 n - 9co 2 ) (4.66) 
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C 3n = n I ” (3/2) coKI® n + 3coKI« n - ( 5/2) Al^/g 2 

+ (i/4)wKI^/g 2 - (1/16)K 8 I 3 “ - ( i/16) K 3 Ig^ 

- (5/16)K 3 co 8 I 3 ^/g 4 - (1/16) K 3 Ig 4 - (3/8)o> 4 K 3 I^/g 2 
+ (i/8)K 3 I^ + (3/8)K 3 co 4 I^/g 2 ]/(a;J n - 9c o 2 ) 

(4.67) 


All values of 


L 

c 

qn 


equal zero except the ones given by equations (4. 64) 


through ( 4 . 67 ) . 


d 1 = 0 (4.68) 

qn 

B 2 = TTlwKl/i + (i/2)wKI 1 5 1 + wKlJ - ( 1/2) c^KI^/g 2 
-co 3 K if/g 2 + ( 1/2) co K ii/g 2 + ( 1/4) w K I^/g 2 
+ (9/16)K 3 lt} - (3/t6)K 3 lu - (3/16) K 3 co 8 I if/ g 4 

- ( 3/ 1 6) K 3 lH + (3/8) K s co 4 liVg 2 + (3/8) K 3 1}} 

+ (l/8)K 3 o; 4 I 1 1 1 7 /g 2 ]/(2a;K) (4.69) 

D 2 = 0 (4.70) 


The third-order term in the expansion of the velocity potential can 
now be expressed as 


= E l< c ln ,t ’in + c L^J cos(t > 


n=i 


3n r 3n' 


+ (C m^ln + C 3n^3n )COS(3T)] 


(4. 71) 
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The wave height rj 3 can be evaluated from equation (2.25d) once equations 
(4.15), (4.16), (4.17), (4.25), (4.45), and (4. 71) have been used. 

A nonlinear solution has been developed through the third-order term, 
but two unknown quantities, Kande, still exist. To evaluate K and obtain a 
meaningful physical significance for e , consider the first approximation to 
the wave height rj . 

7] = erj 1 = eKoj sin(T) 4>u/g on S 0 (4.72) 

Since K is arbitrary, let it be evaluated so that 

r? 1 = 1 at r= ir/2 , 0=0, and r=R ; (4.73) 

therefore K becomes 

K = g/[o> <^(11,0,0)] (4.74) 


and 


r] = e at r = n/2 , 0=0, and r = R ( 4. 75) 

Thus the wave height of the linear theory evaluated atT= 7 r/ 2 , 0=0, and 
r = R is equal to e . This is the same physical significance found by Mack 
[4] , DiMaggio and Rehm [5] , and Baird [6] . 

It is useful to manipulate the nonlinear frequency f so that a 
frequency correction factor can be extracted from it. This is accomplished 
by squaring f and rearranging the terms 
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(t ) 2 = w 2 [i + ( 2B 2 / co ) e 2 + . . . ] 


(4. 76) 


Thus, for a solution through third-order terms, the nonlinear frequency 

parameter (r) 2 is different from the linear frequency parameter by a factor 

of 2B 2 /co , which is called the frequency correction factor G . It is also 

c 

convenient to let 

(t) 2 = X* (4.77) 

and 

co 2 = A (4.78) 

Then equation ( 4. 76) , for a third-order solution, becomes 

A* = A ( 1 + G e 2 ) , (4.79) 

c 

where 

G = 2B 2 /o> (4.80) 

c 

The solution of the nonlinear fluid oscillation problem in an axisym- 
metric container has been formulated, but some things of interest should be 
noted. In the expansion of K , all the terms were zero. One would expect 
this since K is an amplitude parameter and should remain constant for a con 
servative system such as this one. In the expansion of r, it was found that 
B 1 was zero. Thus, the second -order equations do not affect the frequency 
of oscillation. 
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NUMERICAL RESULTS 


Numerical results were obtained for finite -amplitude standing waves 
near the first antisymmetric linear mode in three axisymmetric containers. 
Since cylindrical tanks with bulkheads of various shapes are in general use, 
it was considered appropriate to apply the theory developed in this work to 
obtain numerical results for them. The containers for which results were 
obtained were a cylindrical tank with an ellipsoidal bulkhead, a cylindrical 
tank with a conical bulkhead, and a cylindrical tank with a truncated conical 
bulkhead. 

In the case of the cylindrical tank with an ellipsoidal bulkhead, the 
eccentricity of the bulkhead was considered a parameter and varied accord- 
ingly. The eccentricity, e , is taken to mean the ratio of the depth to the 
radius of the bulkhead, and numerical results were obtained for values of 
e = 0.5, e = i.O, and e = 2.0. 

The parameter that was varied in the study of the cylindrical tank with 
a conical bulkhead was the angle between a horizontal line and the generator 
of the cone. This angle /3 took on values of 30 deg, 45 deg, and 60 deg. 
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The third tank geometry studied is similar to the conical bulkhead 
tank except that the bulkhead is truncated at one -half the depth of the cone. 
Thus the cylindrical tank, with a truncated conical bulkhead, has a radius at 
the bottom equal to one -half the radius of the cylindrical portion of the tank. 
Again, the angle /3 was parameterized for values of 30 deg, 45 deg, and 60 
deg. 

It is convenient to nondimens ionalize the numerical results in such a 
way that one is able to compare the effect of the various tank geometries on 
the results. The two dimensional quantities in this analysis are length and 
time. Thus, the radius of the cylindrical portion of the tanks, R , and the 
local gravitational acceleration, g , will be used for this purpose . The 
nondimensional variables of interest are 

H = H/R , e = e , A = AR^g , A* 

|3 = ]8 , G c = G c R^ , e = e/R t , r? 

Linear Resuits 

To obtain the nonlinear solution it is necessary to solve the linear 
problem first. Thus, the accuracy of the method by which the linear results 
were achieved must be considered. 

It was explained previously that the axisymmetric container being 
analyzed was generated by rotating a curve y about the axis of symmetry, 
and this curve was expressed in parametric form. It was assumed that in 


= A*R t /g , 

= n/R t 
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writing a computer program for this problem, the curve y was made up of 
several sections, and each of the individual curves was of the form 


r = a. 
i 

cos(b.a 

+ C i } 

+ d.a 2 
i 

+ e.a 

i 

+ f. 
1 

(5.1a) 

z = g. 

sin^a 

+ P i> 

+ q.Q! 2 
i 

+ S.Q! 
1 

+ t. 
1 

(5.1b) 


where a is the parameter. 

Thus, by using several sets of equations and selecting the proper 
coefficients in equation (5. 1) , an axisymmetric container of quite a general 
shape can be obtained. 

The linear solution is obtained by solving the eigenvalue problem, 
which is represented by equation (3. 21) and from which a great deal can be 
learned about the convergence of the solution. Symbolically, equation (3.21) 
can be written as 

{[A] + [B] - 2u 2 /g[C] }{a} = 0 (5.2) 

The assumption was made previously that the velocity potential of an 

arbitrary axisymmetric container could be approximated by a series of 

velocity potentials from a cylinder of radius R and height H , where 

c c 

R equals the maximum radius of the arbitrary container, and H equals 
c c 

its height. Thus for the containers under study, R is the radius at the 

c 

free surface. Bearing this in mind, matrices [B] and [C] are diagonal. 
Matrix [A] is obtained by integrating a function of the assumed velocity 
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potentials, which is zero at R , over the container walls. Therefore, if 

c 

the container being studied was a cylinder, this function would be zero at the 

container walls, and [A] would become zero. Thus matrix [A] can be 

called a deviation matrix since it is a measure of the degree to which the 

container deviates from a cylinder. A careful study of the function that is 

integrated to obtain [A] can explain something about the convergence of the 

problem. The function is zero at R and becomes larger as r is decreased; 

c 

also, the function decreases with fluid depth. Thus one would expect excellent 

convergence for a container whose radius deviated only slightly from R near 

c 

the free surface, even if it was greatly different at larger depths. 

A study was made varying the number of terms, NT , in the approxi- 
mating series as a check of the convergence. The containers used in this 
study were a cylindrical tank with an ellipsoidal, a conical, and a truncated 
conical bulkhead. Each was filled to a depth of H = 1. 0 . Tables I, II, and 
III show the resulting eigenvalues to six decimal places. As expected, the 
best convergence was observed in the cylindrical tank with the truncated 
conical bulkhead since it was filled above the conical section, and the devia- 
tion matrix should have been small. The worst convergence was observed 
for the conical bulkhead tank since the fluid was in the conical section. 

Figure 2 shows a comparison of the first three nonaxisymmetric 
frequencies of a fluid in a cylindrical tank with a spherical bulkhead with 
computed and test results obtained by Budiansky [13] . The comparison is 

excellent, indicating that the method by which the linear solution was obtained 
is valid. 
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Table I. Linear Frequencies of a Fluid Contained in a Cylindrical Tank 
with an Ellipsoidal Bulkhead for Nine Values of NT 



*k2 

1.57883 

5.28157 

1.57868 

5.28119 

1.57865 

5.28108 

1.57863 

5.28104 

1.57863 

5.28102 

1.57862 

5.28101 

1.57862 

5.28101 

1.57862 

5.28101 

1.57862 

5.28100 


\8 Pk9 


1.57862 5.28101 8.50646 11.6845 14.8468 18.0018 21.1527 24.3013 


1.57862 5.28101 8.50646 11.6845 14.8468 18.0018 21.1527 24.3012 27.4482 


1.57862 5.28100 8.50646 11.6845 14.8468 18.0017 21.1527 24.3012 27.4481 30.5940 














































































Table II. Linear Frequencies of a Fluid Contained in a Cylindrical Tank 
with a Conical Bulkhead for Nine Values of NT 


— — — — — — — — — — — — — — — “ — — — 

k = 1 H = 1 0 = 45° 

NT 

\l 

\2 

\.S 

\4 

\5 


\7 

\s 


\i0 

2 

1.07943 

4. 85539 









3 

1.06931 

4.79928 

8.10620 








4 

1 . 06444 

4.77214 

8. 10531 

11.2961 







5 

1.06168 

4.75637 

8.02401 

11.2460 

14.4670 






6 

1.05989 

4.74610 

8.00515 

11.2169 

14.4191 

17.6286 





7 

1.05866 

4.73897 

7. 99198 

11.1971 

14.3906 

17.5826 

20.7850 




8 

1.0577 

4. 73377 

7.98230 

11.1825 

14.3706 

17.5549 

20. 7406 

23. 9383 



9 

1.05710 

4.72984 

7.97494 

11.1715 

14.3555 

17.5351 

20.7137 

23.8952 

27.0894 


10 

_ .. J 

1.05657 

4. 72676 

7. 96916 

11.1628 

14.3438 

17.5199 

20.6942 

23. 8690 

27. 0474 

30.2389 




g Table III. Linear Frequencies of a Fluid Contained in a Cylindrical Tank 

with a Truncated Conical Bulkhead for Nine Values of NT 


k = l H = i j8 = 45° 

NT 


*k2 

*k3 

\4 

\l5 

\:6 

*k7 

^k8 

\9 

\l0 

2 

1.68639 

5.32880 

' 








3 

1.68639 

5.32880 

8.53622 








4 

1.68639 

5.32880 

8.53622 

11.7060 







5 

1.68639 

5.32880 

8.53622 

11.7060 

14. 8636 






6 

1.68639 

5.32880 

8.53622 

11.7060 

14. 8636 

18.0155 





7 

1.68639 

5.32880 

8.53622 

11.7060 

14. 8636 

18.0155 

21.1644 




8 

1.68639 

5.32880 

8. 53622 

11.7060 

14. 8636 

18.0155 

21.1644 

24.3113 



9 

1.68639 

5.32880 

8.53622 

11.7060 

14. 8636 

18.0155 

21.1644 

24.3113 

27.4570 


10 

1.68639 

5.32880 

8.53622 

11.7060 

14. 8636 

18.0155 

21. 1644 

24.3113 

27.4570 

30.6019 




LINEAR FREQUENCY 



FLUID DEPTH - H 


Figure 2. Comparison of Linear Frequencies in a Cylindrical Tank 
with a Spherical Bulkhead as a Function of Fluid Depth 
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Figure 3 shows the linear frequency of a fluid in a cylindrical tank 
with an ellipsoidal bulkhead as a function of depth. Curves are presented for 
e = 0.5, e = 1 . 0 , and e = 2. 0 . 

Figure 4 shows the linear frequency of a fluid in a cylindrical tank 
with a conical bulkhead as a function of depth. Curves are presented for 
|3 = 30 deg , /3 = 45 deg , and n = 60 deg . As can be observed, the 
frequency decreases with depth until the fluid level reaches the conical 
section, and then it starts to increase. 

Figure 5 shows the linear frequency of a fluid in a cylindrical tank 
with a truncated conical bulkhead as a function of depth. Curves are 
presented for /3 = 30 deg , j3 = 45 deg , and /3 =60 deg . In comparison 
with Figure 4, one notes that they are similar until shallow depths are 
reached. This difference, of course, is because the flat bottom of the 
container becomes effective. 

It should be noted that as the fluid depth increases, the tank bottom 
becomes less effective. Thus, as the depth increases, the frequency of the 
fluid in all three tanks approaches the frequency obtained for a fluid in a 
cylindrical tank of infinite depth. 

Nonlinear Results 

The assumption, w?. ^ i 2 co 2 , was made previously. At certain 
critical depths this assumption is violated, and the solution becomes invalid. 
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LINDAR FREQUENCY 



0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 1.6 1.8 2.0 

FLUID DEPTH - H 


Figure 3. Linear Frequency of a Fluid in a Cylindrical Tank 
with an Ellipsoidal Bulkhead as Function of Depth 
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LINEAR FREQUENCY 



0.0 0.2 s 0.4 0.6 0.8 1.0 1.2 1.4 1.6 1.8 2.0 

FLUID DEPTH -H 


Figure 4. Linear Frequency of a Fluid in a Cylindrical Tank 
with a Conical Bulkhead as a Function of Depth 
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LINEAR FREQUENCY 






Mack [ 4] has pointed out the reason for the solution becoming invalid at the 
critical depths. It has been assumed that there is a first mode of order e 
oscillating at a frequency co and that all other modes and harmonics are of 
order e 2 or higher. However, at the critical depths this assumption is 
invalid. In general, critical depths do not cause any problem since the 
solution is known to be invalid and is disregarded. However, at shallow 
depths it is impossible to do this. Thus, for all cases studied, the 
solutions were found to be invalid at shallow depths. Therefore, an 
experimental study should be performed to determine the range for 
which the solution is valid. 

In an attempt to verify this analysis, the frequency correction factor 

G was computed for a cylindrical tank as a function of depth and compared 
c 

with the solution obtained by DiMaggio and Rehm [ 5] as shown in Figure 6. 
As can be observed, the comparison is excellent. 

Figure 7 shows a comparison of the nonlinear and linear frequency 
of a fluid in a cylindrical tank with a spherical bulkhead as a function of 
depth. It can be observed that the nonlinear frequency is lower than the 
linear frequency above a depth of H = 0. 35 and higher below this depth. 

The linear and nonlinear wave profile was computed for e = 0. 2 and 
presented in Figure 8. The curve is shown in a configuration of maximum 
potential energy, which occurs at t = [ (1/2) + n] 7 r and which is, of course, 
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FREQUENCY CORRECTION FACTOR - 



Figure 6. Comparison of Frequency Correction Factor 
for a Cylindrical Tank as a Function of Depth 
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FREQUENCY 



Figure 7. Comparison of Nonlinear and Linear Frequency 
of a Fluid in a Cylindrical Tank with a 
Spherical Bulkhead as a Function of Depth 
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VfAVE HEIGHT 



1.0 0.8 0.6 0.4 0.2 0.0 0.2 0.4 0.6 0.8 1.0 

RADIUS - r 


Figure 8. Comparison of Linear and Nonlinear Wave Heights 
in a Cylindrical Tank with a Spherical 
Bulkhead as a Function of Radius 
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when the velocity is zero. The nonlinear wave profile is different from the 

linear wave in that a nodal line does not exist and the surface is never flat. 

Figure 9 shows the frequency correction factor for a fluid in a 

cylindrical tank with an ellipsoidal bulkhead as a function of depth. It is 

interesting to note that for e = 0. 5 and e = 1. 0 the curve for G is similar 

c 

in shape to the one for a cylinder, but for e = 2. 0 it is quite different. 

Observing Figure 3, one notes that for e = 2. 0 the linear frequency is 

increasing with a decrease in depth until the fluid level becomes small. It 

appears that an increase in linear frequency with a decrease in depth tends to 

make G decrease . 
c 

Figure 10 shows the frequency correction factor for a fluid in a 

cylindrical tank with a conical bulkhead as a function of depth. For all three 

values of /3 , G c tends to increase with a decrease in depth until the fluid 

surface is in the conical bulkhead, at which time it starts to decrease . This 

can be explained as before by the fact that the linear frequency is increasing 

(Fig. 4) as the depth is decreasing. 

Figure 11 shows the frequency correction factor for a fluid in a 

cylindrical tank with a truncated conical bulkhead as a function of depth. 

For /3 = 30 deg , G behaves very much like it does for a cylindrical tank, 

c 

which is expected since the solution becomes invalid before it becomes small 

enough for the fluid to be in the bulkhead section of the tank. It is seen that 

G for (3 = 45 deg and /3 = 60 deg tend to increase with a decrease in 
c 
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0,0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 1.6 1.8 2.0 

FLUID DEPTH . H 


Figure 9. Frequency Correction Factor for a Fluid in a 
Cylindrical Tank with an Ellipsoidal Bulkhead as a Function of Depth 
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CONCLUSIONS 


A solution for finite -amplitude free oscillations of a fluid in a partially 
filled axisymmetric container has been presented. The formulation of the 
problem resulted in a nonlinear boundary value problem where the nonlinearity- 
occurred in the boundary condition at the free surface of the fluid. The 
boundary condition at the container wall and the differential equation were 
linear. 

The solution was obtained by first linearizing the free surface boundary 
condition and solving the resulting linear boundary value problem. Then the 
linear solutions, which satisfy the differential equation and the boundary 
condition at the container wall, are used to construct a solution that satisfies 
the nonlinear boundary condition at the free surface asymptotically. A solu- 
tion through the third-order term was developed. 

Numerical results were obtained for finite -amplitude standing waves 
near the first antisymmetric linear mode in three axisymmetric vessels. 
Results were found for a cylindrical tank with an ellipsoidal bulkhead, a 
cylindrical tank with a conical bulkhead, and a cylindrical tank with a 
truncated conical bulkhead. 
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The method of solving the linear problem was verified by comparing 
the computed linear frequency with published test results for the case of a 
cylindrical tank with a hemispherical bulkhead. The nonlinear solution was 
checked by comparing the results for a cylindrical tank with published 
theoretical results. 

George C. Marshall Space Flight Center 

National Aeronautics and Space Administration 

Marshall Space Flight Center, Alabama 35812, August 29, 1969 
981-10-10-0000 
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